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Abstract. We have used the available samples of Gigahertz Peaked Spectrum (GPS) galaxies to investigate their 
evolution properties in the framework of the "youth" scenario. Care was taken to properly allow for the different 
selection criteria used to define the samples. We find that the observed redshift and peak frequency distributions 
can be satisfactorily accounted for in terms of simple luminosity evolution of individual sources, along the lines 
discussed by Fanti et al. (1995) and Begelman (1996, 1999), although the derived parameter values have large 
uncertainties due to ambiguities in the selection of GPS sources and to the incompleteness of redshift measure- 
ments. However the simplest self-similar model, whereby the evolution is controlled only by the radial profile of 
the density of the ambient medium is not good enough and one additional parameter needs to be introduced. The 
fit requires a decrease of the emitted power and of the peak luminosity with source age or with decreasing peak 
frequency, at variance with the Snellen et al. (2000) model. On the other hand, our analysis confirms the rather 
flat slope of the luminosity function, found by Snellen et al. (2000) who also report indications of a high luminosity 
break, not required by the data sets we have used. Our results suggest that the GPS galaxies are the precursors 
of extended radio sources with luminosities below the break of the luminosity function. No cosmological evolution 
of the GPS galaxy population is required by presently available data. 
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1. Introduction 

Several lines of evidence increasingly point to Gigahertz 
Peaked Spectrum (GPS) sources (see O'Dea 1998 for a 
comprehensive review) being extremely young. Spectral 
age constraints, energy supply arguments and, most con- 
vincingly, VLBI kinematic studies, are all consistent with 
the current radio activity in these objects having turned 
on less than about 1000 years ago (Conway 2002), al- 
though in some cases there is evidence that this is not 
the first instance of activity, but sources appear to have 
been re-activated after a period of quiescence (Baum et al. 
1990), perhaps by effect of merger activity (Tingay et al. 
2003). These sources thus offer the exciting opportunity 
of studying radio sources in early stages of development. 

There is a clear anti-correlation between the peak 
(turnover) frequency, v p , and the projected linear size of 
GPS sources (Fanti et al. 1990; O'Dea & Baum 1997). 
Although this anti-correlation does not necessarily define 
the evolutionary track, a decrease of the peak frequency as 
the emitting blob expands is indicated, suggesting a sce- 
nario whereby GPS sources evolve into Compact Steep 
Spectrum (CSS) sources and, finally, into the classical 
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large scale radio sources (Fanti et al. 1995; Readhead et 
al. 1996; Begelman 1996; Snellen et al. 2000). 

The self-similar evolution model by Begelman (1996, 
1999) imply that the evolution of the radio power as the 
source expands, depends on the radial profile of the exter- 
nal density. To test this scenario and to constrain the birth 
rates and the evolutionary properties of these sources, de- 
tailed comparisons of the model predictions with survey 
data are necessary. An exploratory study in this direction 
was presented by De Zotti et al. (2000), who showed that 
this kind of models may account for the observed counts, 
redshift and peak frequency distributions of the samples 
then available, but for rather unexpected values of the 
parameters. 

Since then, new important surveys of GPS sources 
have been published, extending the coverage of the peak 
frequency-redshift plane. The Dallacasa et al. (2000) 
and Bolton et al. (2004) samples comprise sources with 
high turnover frequencies (called High Frequency Peakers, 
HFPs), which are presumably the youngest. Snellen et al. 
(2002) and Edwards & Tingay (2004) selected samples 
of southern GPS sources, thus improving the sky cover- 
age. Snellen et al. (2004) focussed on the relatively rare 
low-z sources. The first year WMAP survey (Bennett et 
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al. 2003b) has provided a complete sample of over 200 
sources with simultaneous flux density measurements at 4 
frequencies (22.8, 33.0, 40.7, 60.8, and 93.5 GHz), that can 
be exploited to search for extreme GPS sources, peaking 
at mm wavelengths. 

Moreover, the variability study of the Dallacasa et 
al. (2000) sample, carried out by Tinti et al. (2005), has 
shown that most HFP candidates classified as quasars are 
most likely blazars caught by during a flare of a highly 
self-absorbed component dominating the emission, while 
candidates classified as galaxies are consistent with being 
bona fide HFPs. Polarization measurements (Dallacasa et 
al., in preparation) lend further support to this conclu- 
sion. Torniainen et al. (2005) from a study of the long 
term variability of 33 objects previously classified as GPS 
or HFP sources, mostly identified with quasars and mostly 
peaking above 5 GHz, concluded that only 5 keep the GPS 
properties over time. Clearly, a serious blazar contamina- 
tion of GPS samples may lead astray analyses of their 
statistical properties. 

In view of these new data, the previous conclusions on 
the evolution of the GPS population need to be recon- 
sidered in depth. In this paper we present a new analysis, 
still in the framework of Begclman's (1996, 1999). Since, as 
mentioned above, there are strong evidences that samples 
of GPS sources identified as quasars are heavily contam- 
inated by flaring blazars, we confine our analysis to GPS 
galaxies. Our reference evolutionary scenario is briefly de- 
scribed in Sect. 2, in Sect. 3 we give an account of the 
data sets we use, in Sect. 4 we analyze the effect of the 
different selection criteria adopted to define the samples, 
in Sect. 5 we present our results, in Sect. 6 we summarize 
and discuss our main conclusions. 

We have adopted a flat ACDM cosmology with f^A = 

0. 73 and H a = 70^8-* Mpc -1 (Bennett et al. 2003a). 

2. The evolutionary scenario 

Following De Zotti et al. (2000), we make the following 
assumptions: 

1. The initial radio luminosity function (in units of 
Mpc~ 3 dlogii -1 ) is described by a power law: 



nil, [ j- 



(i) 



where Li is the luminosity before self-absorption. 
2. In the GPS phase, the properties of the sources are de- 
termined by the interaction of a compact, jet-driven, 
overpressured, non thermal radio lobe with a dense 
interstellar medium (Bcgclman 1996, 1999; Bicknell 
et al. 1997); the timescale of the interaction is very 
short in comparison with the cosmological-expansion 
timescale, so that the luminosity evolution of individ- 
ual sources occurs at constant z. As the radio lobe ex- 
pands in the surrounding medium, the emitted radio 
power varies with the source age, r, as Li oc r~ r ', and 
its linear size I varies as I oc r e . If the density of the 



surrounding medium scales with radius as p e oc r ™, 
we have (Begelman 1996, 1999) r/ = (n + 4)/[4(5 - n)] 
and e = 3/(5 — n). The Begclman's dependence of 
the radio power on source age assumes that the syn- 
chrotron emission is dominated by the high pressure re- 
gion close to the head of the jet. Alternatively, Snellen 
et al. (2000) and Alexander (2000) assume that the 
emission comes from the full volume of the expanding 
cocoon, so that 77 = (8 — 7n)/4(5 — n). In Begelman's 
model the luminosity decreases with increasing source 
size for n > 0. On the other hand, according to Snellen 
et al. (2000) and Alexander (2000) the radio luminos- 
ity increases as l 2 / 3, if n = (i.e. within the galaxy 
core radius) and then decreases as Z -0 - 5 as the density 
profile steepens to n — 2 at larger radii. Evolutionary 
tracks where the luminosity increases within the core 
radius and then decreases in a falling atmosphere in the 
power-diameter plane are also presented by Carvalho 
& O'Dea (2003) based on the results of a series of 
hydrodynamical simulations. There is a clear anti- 
correlation between intrinsic turnover frequency, v p , 
and linear size (O'Dea & Baum 1997): v p oc l~ 5 , 
with 6 ~ 0.65. It follows that v p scales with time as 
v p oc t~ a , with A = 8e. 
3. The spectra of GPS sources are described by: 



L v — Lp x 



(v/v p ) aa if v <v p 
(vjvp)~ a if v > Vp 



(2) 



with a a = 0.8 and a = 0.75, the mean values found by 
Snellen et al. (1998b). 

As the radio lobe expands, the peak luminosity L p varies 
as the consequence of the variation of both the emitted 
radio power and of v p . Hence: 

— = — (3) 

Vp.% J \ f Pt i J 

with L Pt i(z) oc Li(z) and p = —r\ + a\. 

If the birth rate of GPS sources is constant on 
time scales much shorter than the cosmological-expansion 
timescale, the peak luminosity function per unit dlogi p 
is: 

n(Lp) cx L'Jp. (4) 

Also, since, in this case, the comoving number of sources of 
age r within dr is simply proportional to dr and dr/d^ p oc 
(vp/vp^^^ 1 /^, the epoch dependent luminosity func- 
tion at a given frequency v (Mpc~ 3 dlogL„ -1 GHz -1 ) 
writes: 



n(L u , u p , z) = n { 



L Pt i(L u , Vp) 



-0 



-(l+l/A) 

P ) (5) 



where L±(z) is the redshift-dependent normalization lu- 
minosity. We have assumed luminosity evolution and 
adopted a very simple parameterization for it: 



(1 -|- z\ ifz<Cz 
L *( z ) = L ox{ {1 + Zc) k iiz>z c c 



(6) 
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The redshift z c at which luminosity evolution levels off is 
a model parameter. We have normalized monochromatic 
luminosities to L a — 10 32 ergs -1 Hz -1 . 

The luminosity function at a frequency v 2 is related to 
that at a frequency v\ by: 



Table 1. Summary of GPS samples. 



n(L V2 ) = n(L Vl )x ■ 



if v p < v\ < v 2 



The number counts per steradian of GPS sources brighter 
than S u at the frequency v, with an observed peak fre- 
quency max(v, i/p, mi n) < Vpfi < t^max are given by 

rraixL[z f ,z m (S u )] ^y 

N (> S v \ v p<0 > v) = / dz — x 

./n dz 



X 



dup x 



/ max[z/ p . m i n (l+z),i/(l+z)] 
/•log L max (z) 

x / dlogij, n{L v , v p , z) . 

•/log L m j n (S„ 



(8) 



where Zf is the redshift of formation of the first GPS 
sources, z m is the maximum redshift at which sources can 
have a flux > S v , £ m in is the minimum luminosity of a 
source of given z and v p yielding a flux > S v , dV/dz is 
the comoving volume element within a solid angle to: 
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dz ~ H U \l + z) 2 E{z) 
cIl(z) being the luminosity distance: 



, / x c , , [ z dz' 



e{z) = (n M (i + zf + n A ) 



1/2 



(9) 



(10) 



(11) 



The flux density at the frequency v is related to the rest- 
frame luminosity L v by: 



S v 



L v K{z) 



K(z) = (l + z) 



(12) 
(13) 



K(z) being the K-correction. 

Similarly, the number counts of GPS sources with an 
observed peak frequency Vp^m\ n 
are given by 

r-vain[z f ,z m {S v )] 

N (> S v ] u pfl < v) = / dz — x 

Jn dz 



X 



vci\n[v (l+z),i/ p max (l+z),i/ Pii ] 



dv p x 



„„(l+z) 
/■log L max (z) 

x / d log L v n{L v ,v p , z) , 

•/log L m i n (S„ ,Z,I/ p ) 



(14) 



(yi/vi) Qa if v x < v 2 < v p 

{v\lvp) aa (yilv v ) a if vi < v p < v 2 (7) 



Sample 


•Shin 


V 


area 


V 


max 
P 


tot 




Jy 


GHz 


deg 2 


GHz 


GHz 




HFP 


0.3 


4.9 


15840 


5.34 


11.1 


5 


Stan 


1.0 


5.0 


24600 


0.4 


4.9 


19 


B_SA 


0.025 


15.0 


176 


3.43 


127.0 


10 


B_SB 


0.060 


15.0 


70 


3.29 


10.3 





Snel 


0.018 


0.325 


522 


1.0 


5.7 


14 


CORA 


0.1 


1.4 


2850 


0.460 


2.3 


6 


Park 


0.5 


2.7 


12802 


0.4 


5.0 


48 


Edwa 


0.95 


5.0 


16414 


0.7 


7.5 


9 



The distribution of observed peak frequencies per unit 
dvpfl in a flux limited sample, M(v P fi\ > Su), is given by: 

rIDm[Zf,Z m (Sv,V p ,o)] fly 

Af (v p . ;> S v )= dz ~r~ x 

Jo dz 

/•log L max (z,v Pi0 (1+z)) 

x / dlogL„n[L„,z,v pfi (l + z)]. (15) 

Jlog L min (S,,,z,u P:0 (l+z)) 



3. GPS samples 

There is no clear-cut definition of GPS sources and differ- 
ent criteria have been adopted for identifying them. On 
the other hand, we need to combine at least the most 
reliable samples to gather enough data for a meaningful 
analysis to be possible. We now describe the samples used 
and the corrections applied to make them as homogeneous 
as possible. The main properties of the samples are sum- 
marized in Table HI where HFP, Stan, B_SA, B.SB, Snel, 
CORA, Park, and Edwa denote, respectively, the samples 
by Dallacasa et al. (2000, HFP sample), by Stanghellini et 
al. (1998), by Bolton et al. (2004, 9C samples A and B), 
by Snellen et al. (1998, as revised by Snellen et al. 2000, 
faint WENSS sample), by Snellen et al. (2004, CORALZ 
sample), by Snellen et al. (2002, Parkes sample), and by 
Edwards & Tingay (2004, ATCA sample). 

3.1. The 9C samples 

Bolton et al. (2004) selected two complete samples of 
sources from the first three regions of the 9C survey: a 
deeper sample (sample A), complete to 25 mJy and con- 
taining 124 sources over a total area of 176 square degrees; 
a shallower sample (sample B), complete to 60 mJy, com- 
prising 70 sources in an area of 246 square degrees, includ- 
ing the area covered by sample A. Restricting the sample B 
to the 70 sq. deg. not overlapping the sample A we are left 
with 31 sources. Simultaneous observations of each source 
were made at frequencies of 1.4, 4.8, 22 and 43 GHz with 
the VLA and at 15 GHz with the Ryle Telescope. In addi- 
tion, 51 sources were observed within a few months at 31 
GHz with the Owens Valley Radio Observatory (OVRO) 
40m telescope. Sources with spectral index between 1.4 
and 4.8 GHz af f < —0.1 (S u oc v~ a ) were referred to as 
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Table 2. GPS candidates from Bolton et al. 2004. 



Name 






Sp 

(mJy) 


Up 

(GHz) 


J0003 + 2740 A 


Q? 




76 ± 1 


5.4 ±0.1 


J0003 + 3010 A 


G? 


0.55 


60 ±3 


10.3 ±0.3 


J0010 + 2854 A 


Q? 




100 ±5 


52 ±6 


J0012 + 3353 A 


G? 


0.50 


199 ± 12 


127 ± 36 


J0012 + 3053 A 


G? 


0.45 


27 ± 1 


20 ±4 


J0020 + 3152 A 


- 


1.1 


43.8 ±0.8 


4.91 ±0.09 


J0024 + 2911 A 


Q? 




42 ±2 


13.9 ±0.5 


J0032 + 2758 A 


G 


0.51 


34.4 ± 0.4 


4.6 ±0.2 


J0919 + 3324 s 


G? 


0.35 


430 ± 16 


10.3 ±0.5 


J0925 + 3127 s 


G? 


0.26 


135 ± 2 


3.29 ±0.06 


J0931 + 2750 s 


G? 


0.49 


169 ±3 


8.3 ±0.3 


J0935 + 2917 A 


Q? 




47.3 ±0.9 


5.7 ±0.2 


J0936 + 3207 A 


G 


0.20 


56 ±5 


14 ± 1 


J0936 + 2624 s 


Q? 




179 ± 18 


3.8 ±0.3 


J0940 + 2603 s 


G? 


0.37 


495 ± 10 


6.1 ±0.8 


J0945 + 3534 s 


Q? 




413 ±8 


4.0 ±0.2 


J0952 + 3512 s 


Q? 




427 ±8 


4.1 ±0.2 


J0955 + 3335 s 


Q? 




107 ±2 


5.9 ±0.2 


J1506 + 4239 A 


G 


0.38 


766 ± 19 


14 ± 1 


J1517 + 3936 A 


G 


0.46 


43 ±2 


21 ±3 


J1521 +4336 A 


Q? 




431 ± 10 


6.2 ±0.2 


J1526 + 3712 A 


G? 


0.61 


74 ±2 


6.9 ±0.3 


J1526 + 4201 A 


G 


0.36 


67 ± 1 


8.1 ±0.2 


J1528 + 3816 A 


G? 


0.52 


80 ±3 


26 ±5 


J1530 + 3758 A 


G 


0.19 


141 ±3 


3.43 ±0.06 


J1540 + 4138 A 


G 


0.17 


46 ±8 


8.8 ±0.2 


J1550 + 4536 A 


G 


0.50 


62 ± 1 


3.5 ±0.07 


J1554 + 4350 A 




1.2 


45.0 ±0.8 


10.9 ±0.3 


J1554 + 4348 A 


G? 


0.36 


63 ± 1 


3.7 ±0.1 


J1556 + 4259 A 


Q? 




94 ±4 


4.3 ±0.2 



GPS sources. The adopted criterion implies that sources 
peaking at v p ^ 5 GHz were preferentially (but not exclu- 
sively) selected. With this definition, there are 22 GPS 
sources in sample A (14 galaxies, 6 quasars and 2 uniden- 
tified sources) and 8 (4 galaxies and 4 quasars) in the 
(redefined) sample B. 

In order to estimate the peak flux densities, S p , and 
frequencies, v p , of the sources, we have fitted the radio 
spectra with the hyperbolic function used by Tinti et al. 
(2005): 

log(S) = log(S p ) + b - [b 2 + c 2 (\ogv- log v p ) 2 } 1/2 . (16) 

The best fit values of v p and S p , obtained minimizing the 
chi-square function with the Minuit package (CERN li- 
braries), are reported with their errors in column 4 and 5 
of Table H 

No redshift measurements are available for these 
sources. However, as shown by Snellen et al. (1996, 2002), 
GPS galaxies show a well defined i?-band Hubble diagram, 
with a low dispersion. In terms of Gunn r magnitudes, 
Snellen et al. (1996) found a best fit relation: 

r = 22.7 + 7.4 log(z). (17) 



Table 3. The HFP galaxy sample. The * denotes an 
empty field, which was attributed z = 1.5. 



Name 


z 


Up 






GHz 


0428±3259 


0.3 


7.3 


0655±4100 


0.02156 


7.8 


1407±2827 


0.0769 


5.34 


1511+0518 


0.084 


11.1 


1735+5049 


1.5* 


6.4 



The conversion from the Kron-Cousins R c magnitudes 
measured by Bolton et al. (2004) to the Gunn r mag- 
nitudes was made taking r = R c + 0.3 (Fukugita et al. 
1995). The redshift estimates are given in Table We 
do not count as galaxies the optically point-like objects, 
classified as G? by Bolton et al. (2004) because of their 
red colours (O — R> 1.6): if they were bright galaxies at 
the redshifts estimated from Eq. (|17f) they should be re- 
solved. This leaves 10 GPS galaxies in sample A and zero 
in sample B (see Table H). 

It is important to note that other criteria for select- 
ing GPS samples include tighter constraints on the low- 
and/or high-frequency spectral indices. The effect of the 
different selection criteria is discussed in Sect. SI 

3.2. The HFP sample 

The bright HFP sample by Dallacasa et al. (2000) was se- 
lected by cross-correlating the 87GB (Gregory et al. 1996) 
sources with S4.9GHZ > 300 mJy with the NVSS catalogue 
(Condon et al. 1998) at 1.4 GHz and picking out those 
with inverted spectra (a < —0.5, S oc v~ a ). It was then 
"cleaned" by means of simultaneous multifrequency VLA 
observations, leaving 55 sources whose single-epoch radio 
spectrum peaks at frequencies ranging from a few GHz to 
about 22 GHz. The sample of HFP candidates comprises 
11 galaxies (including a type 1 Seyfert), 36 quasars, and 
8 still unidentified sources (Dallacasa et al. 2002), over an 
area of 15,840deg 2 . 

Although there are, in this sample, some sources with 
v p < 4.9 GHz, the selection criterion biases the sample 
against such values of v p in a way that we are unable to 
quantify. Therefore we have chosen to confine ourselves to 
sources with v p > 4.9 GHz. Moreover we have excluded 
from the sample the objects with 4.9 GHz flux densities 
smaller than the completeness limit of 300 mJy when were 
they re-observed by Dallacasa et al. (2000) and Tinti et al. 
(2005). After having applied these additional constraints, 
we are left with 5 HFP galaxies, listed in Table |21 

3.3. The faint GPS sample from WENSS 

The selection of this sample is described in detail in 
Snellen et al. (1998). Snellen et al. (2000) applied stricter 
criteria allowing a better control of selection effects; they 
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Table 4. The faint GPS galaxies from WENSS. The * de- 
notes empty fields, which were attributed z = 1.5. Values 
of z with 3 significant digits are spectroscopic, the others 
are photometric estimates. 



Name 


Z 


Vp 








J3U4UU+0U4z 


1 C + 

1.0 


J..U 


T3C\AOd 1 d~\ CO 

J3(J4<3D+OlOZ 


1 c + 
1.0 


l.U 


-dUOoO+D ( 4o 


1 . 


r >j 

. 1 


T2nc;Qn i conn 


1.4 


1.9 


JJUOOUT'-'OIO 


u.uyo 


1 fi 

-L .U 


B1525+6801 


i.i 


1.8 


B1551+6822 


1.3 


1.5 


B1557+6220 


0.9 


2.3 


B1600+7131 


1.5* 


1.7 


B1622+6630 


0.201 


4.0 


B1655+6446 


1.5* 


1.0 


B1841+6715 


0.486 


2.1 


B1942+7214 


1.1 


1.4 


B1946+7048 


0.101 


1.8 



Table 5. The bright GPS galaxies from Stanghellini et al. 
(1998). The * denotes photometric redshift estimates. 



Name 


2 


u p 






GHz 


0019-000 


0.305 


0.8 


0108+388 


0.669 


3.9 


0316+161 


1.2* 


0.8 


0428+205 


0.219 


1.0 


0500+019 


0.583 


2.0 


0710+439 


0.518 


1.9 


0941-080 


0.228 


0.5 


1031+567 


0.459 


1.3 


1117+146 


0.362 


0.5 


1323+321 


0.369 


0.5 


1345+125 


0.122 


0.6 


1358+624 


0.431 


0.5 


1404+286 


0.077 


4.9 


1600+335 


1.1* 


2.6 


1607+268 


0.473 


1.0 


2008-068 


0.7* 


1.3 


2128+048 


0.99 


0.8 


2210+016 


1.0* 


0.4 


2352+495 


0.237 


0.7 



kept only the 14 objects with inverted spectra between 325 
and 5 GHz, and with 325 MHz flux densities > 20 mJy, 
over an area of 522 deg 2 . The redshift of identified sources, 
all classified as galaxies, had to be estimated from their 
optical magnitudes, using Eq. I|17l) . For the 4 unidentified 
sources, also assumed to be galaxies, a redshift of z=1.5 
was assumed. The relevant data for all the 14 objects are 
given in Table 01 



Table 6. The ATCA sample of GPS galaxies from 
(Edwards & Tingay 2004). The * denotes empty fields, 
which were attributed z = 1.5. 



Name 


z 


P 






GHz 


J0241-0815 


0.004 


7.5 


J1543-0757 


0.172 


0.7 


J1658-0739 


1.5* 


4.8 


J1726-6427 


1.5* 


1.1 


J1723-6500 


0.014 


2.7 


J1744-5144 


1.5* 


1.0 


J1939-6342 


0.183 


1.4 


J2257-3627 


0.006 


2.7 


J2336-5236 


1.5* 


1.1 



3.4. The bright GPS sample 

Stanghellini ct al. (1998) selected candidate radio bright 
GPS sources from the Kuhr et al. (1981) catalogue 
(S5GHZ > 1 Jy), over an area of about 24, 600 deg 2 . The 
sample was then cleaned by means of multifrcqucncy VLA 
and WSRT observations, supplemented with literature 
data. They picked out GPS candidates with a turnover 
frequency between 0.4 and 6 GHz, and an optically thin 
spectral index athin > 0.5 (S u tx v~ a ), beyond the peak. 
The final complete sample consists of 33 GPS sources, 
19 of which are identified with galaxies. Four galaxies do 
not have spectroscopic redshift; estimates by Snellen et 
al. (2000) from their optical magnitudes [Eq. l(T7|l ] are de- 
noted by a * in Table [SI where the relevant data for the 
19 galaxies are listed. 

3.5. The ATCA GPS sample 

Edwards & Tingay (2004) have used data from an 
Australia Telescope Compact Array (ATCA) program of 
multi-frequency, multi-epoch monitoring of the portion of 
the VSOP survey sample (Hirabayashi et al. 2000) with 
declinations < 10°. The original sample is defined by: 
S5GHZ > 0.95 Jy a < 0.45, \b\ > 10°. Taking into account 
the further constraint S < 10°, we estimate that the area 
covered is ~ 5 sr. The selected sources have athin > 0.5 
and spectral curvature athin — athick > 0.6. 

We have excluded from the sample the gravitationally 
lensed source J0414 + 0534 because our models do not 
include the effect of lensing, and J1522 — 2730, classified 
as a BL Lac object, whose variability properties suggests 
that it is less likely to be a truly GPS source (Edwards 
& Tingay 2004). This leaves 7 GPS galaxies (5 having 
spectroscopic redshifts), 16 quasars (15 with spectroscopic 
redshift), and 2 empty fields (see Table EJ). 

3.6. The Parkes half- Jy sample 

Snellen et al. (2002) selected from the Parkes multifre- 
quency survey data in a region of about 3.9 sr a south- 
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Table 7. The Parkes half-Jansky sample of GPS galaxies 
(Snellen et al. 2002). The * denotes photometric redshift 
estimates. 



Table 8. The CORALZ sample of local GPS galaxies 
(Snellen et al. 2004). 



Name 


mR 


z 


y 

GHz 


J0022+0014 


18.10 


0.305 


0.6 


J0108-1201 


22.39 


1.0* 


1.0 


J0206-3024 


21.00 


0.65* 


0.5 


J0210+0419 


24.1 


1.5* 


1.3 


J0210-2213 


23.52 


1.4* 


1.5 


J0242-2132 


17.10 


0.314 


1.0 


J0323+0534 


19.20 


0.37* 


0.4 


J0401-2921 


21.0 


0.65* 


1.0 


J0407-3924 


20.40 


0.54* 


0.4 


J0407-2757 


21.14 


0.68* 


1.5 


J0433-0229 


19.10 


0.36* 


0.4 


J044 1-3340 


21.0 


0.65* 


1.2 


J0457-0848 


20.30 


0.52* 


0.4 


J0503+0203 


21.0 


0.583 


2.5 


J0943-0819 


17.50 


0.228 


0.4 


J0913+1454 


20.0 


0.47* 


1.1 


J1044-2712 


21.0 


0.65* 


0.8 


J1057+0012 


21.0 


0.65* 


1.6 


J1109+1043 


20.50 


0.55* 


0.5 


J1110-1858 


19.60 


0.497 


1.0 


J1120+1420 


20.10 


0.362 


0.4 


J1122-2742 


21.0 


0.65* 


0.8 


J1135-0021 


16.50 


0.16* 


0.4 


J1203+0414 


18.80 


0.33* 


0.4 


J1345-3015 


21.0 


0.65* 


2.5 


J1347+1217 


15.20 


0.122 


0.4 


J1350-2204 


20.93 


0.63* 


0.4 


J1352+0232 


20.00 


0.47* 


0.4 


J1352+1107 


21.0 


0.65* 


3.6 


J1447-3409 


21.00 


0.65* 


0.5 


J1506-0919 


19.70 


0.43* 


0.6 


J1521+0430 


22.10 


1.296 


1.0 


J1543-0757 


17.40 


0.172 


0.4 


J1546+0026 


20.10 


0.556 


0.6 


J1548-1213 


21.88 


0.883 


0.4 


J1556-0622 


22.20 


0.94* 


0.4 


J1604-2223 


18.75 


0.141 


0.6 


J1640+1220 


21.36 


1.150 


0.4 


J1648+0242 


21.0 


0.65* 


3.4 


J1734+0926 


20.80 


0.61* 


1.0 


J2011-0644 


21.18 


0.547 


1.4 


J2058+0540 


23.40 


1.381 


0.4 


J2123-0112 


23.30 


1.158 


0.4 


J2130+0502 


22.21 


0.990 


1.0 


J2151+0552 


20.20 


0.740 


5.0 


J2212+0152 


22.0 


0.88* 


0.4 


J2325-0344 


23.50 


1.4* 


1.4 


J2339-0604 


22.91 


1.2* 


0.4 



ern/equatorial sample of GPS sources with 1S2.7GHZ > 
0.5 Jy, excluding objects identified as quasars. The sample 
(see Table 01 consists of 49 objects with spectra peaking 
at v v > 0.4 GHz, 38 of which are identified with galaxies, 
10 are too faint to be identified and 1 is too close to a 



Name 


z 


P 






GHz 


J073328+560541 


0.104 


0.460 


J073934+495438 


0.054 


0.950 


J083139+460800 


0.127 


2.200 


J090615+463618 


0.085 


0.680 


J131739+411545 


0.066 


2.300 


J171854+544148 


0.147 


0.480 



bright star to allow identification, and is excluded from 
the statistical analysis. The authors argue that, based on 
the magnitude distribution of other GPS samples, the 10 
faint sources are unlikely to be quasars, and we assume 
them to be galaxies. Spectroscopic redshifts are available 
for 18 objects. Estimates, or lower limits, for the others 
have been obtained through Eq. i|17|) . No restrictions on 
either the low- or the high-frequency spectral index are 
mentioned. 

3. 7. The CORALZ sample 

The sample of Compact Radio sources at Low Redshift 
(CORALZ) was selected by Snellen et al. (2004) solely 
on the basis of the radio angular size (9 < 2"), indepen- 
dent of radio spectra, picking out sources associated with 
relatively bright galaxies. It is estimated to be ~ 95% 
complete for Si.4GHz > 100 mJy over the redshift range 
0.005 < z < 0.16 in a region of 2850 square degrees. The 
sample comprises 6 GPS galaxies with v v > 0.4 GHz, all 
with spectroscopic redshifts (Table 

3.8. The WMAP survey 

The WMAP point source catalogue (Bennett et al. 2003b) 
comprises 208 point sources with a > 5 cr detection in any 
of the 5 WMAP bands, in the range 22.8-93.5 GHz; 203 
sources have counterparts in existing 4.85 GHz, the re- 
maining 5 sources being probably spurious. 

This is the first simultaneous multifrequency survey 
at mm wavelengths and is therefore well suited to se- 
lect extreme HFPs. Although the completeness limit is at 
522.8GHz ^ 1-25 Jy (De Zotti et al. 2005), to avoid that the 
estimates of spectral indices are too affected by measure- 
ment errors we have confined ourselves to sources with 
S l 22.8GHz > 2Jy, and we have picked out those with in- 
verted spectrum (a < 0) from 22.8 to 33 GHz and from 
4.85 to 22.8 GHz. There are 18 sources satisfying these cri- 
teria, a number close to the prediction of the De Zotti et 
al. (2000) model, for a maximum initial peak frequency of 
200 GHz. Most of these sources are well known calibrators 
and have therefore many observations at many frequencies 
(see Trushkin 2003). They all show strong variability, con- 
sistent with that observed for blazars (see Tinti 2005); 9 of 
them are classified as blazars by Terasranta et al. (1998) or 
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Fig. 1. Peak flux density versus observed peak frequency (left) and peak luminosity versus rest frame turnover 
frequency (right) for GPS galaxies in the various samples. HFP: filled circles; Bolton sample A: x; WENSS: filled 
squares; CORALZ: filled triangles; Stanghellini: *; Parkes: +. 
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Fig. 2. Distributions of measured and estimated red- 
shifts for the HFP, Stanghellini et al. (1998), Snellen et al. 
(1998), Bolton et al. (2004, sample A), Parkes, and ATCA 
GPS galaxy samples. Note that the peaks at z = 1.5 are 
due to unidentified sources, assumed to be galaxies at that 
redshift. 

Donato et al. (2001). The data on the WMAP sample thus 
confirm that most quasars showing a peak at tens of GHz 
are likely blazars caught during a phase when a flaring, 
strongly self-absorbed synchrotron component dominates 
the emission spectrum (Tinti et al. 2005). 

4. Impact of different selection criteria 

The observational properties that generally define the dif- 
ferent samples are: 

— the flux density at the frequency of selection; 



Table 9. Median values of the optically thick and thin 
spectral indices for the various samples. 



Sample 


Qthick 


Qthin 


HFP 


-0.86 


0.61 


Stan 


-0.76 


0.67 


B_SA 


-0.38 


0.76 


B_SB 


-0.43 


0.82 


Snel 


-0.67 


0.64 



— the turnover frequency range; 

— the spectral index in the optically thin and thick parts 
of the spectrum. 

To effectively explore the luminosity and peak frequency 
evolution of GPS sources, we need samples that provide a 
wide coverage of the v p -S p plane. The coverage provided 
by the present samples is shown in Fig. ^ HFP sources 
(Dallacasa et al. 2000) have turnover frequencies greater 
than 4.9 GHz, at or above the upper limits of the samples 
of Stanghellini et al. (1998) and Snellen et al. (1998, 2002, 
2004). Because of their high selection frequency, 15 GHz, 
the Bolton et al. (2004) samples explore the high peak 
frequency region. 

The samples by Bolton et al. (2004) have flux limits 
similar to the faint GPS sample selected by Snellen et al. 
(1998), although at a very different frequency, while the 
HFP and the Stanghellini et al. (1998), and the Parkes 
samples contain bright sources. 

The redshift distributions of the galaxies in the various 
samples are shown in Fig. Many redshifts are estimated 
from optical magnitudes. The peak at z = 1.5 is mostly 
due to optically unidentified GPS sources, assumed, fol- 
lowing Snellen et al. (2000), to have that redshift. 
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Fig. 3. Spectra of GPS sources in the various samples, normalized in both frequency and flux density. 



As noted above, selection criteria include constraints 
on the spectral indices in the optically thick and/or in the 
optically thin spectral region. Stanghellini et al. (1998) 
required athin > 0.5, Dallacasa et al. (2000) demanded 
£*thick < —0.5; Edwards & Tingay (2004) applied both 
constraints. For other samples, the adopted spectral cri- 
teria are less explicit. Snellen et al. (1998, 2000) require 
an inverted spectrum at low frequencies and that the Full 
Width at Half Maximum of the spectrum is less than 2 
decades in frequency; the latter condition implies that, 
typically, athin and —athick are > 0.3, a somewhat less re- 
strictive constraint than adopted for the previously men- 
tioned samples. An even looser criterion was adopted by 
Bolton et al. (2004): af f < -0.1. Not many details are 
given on spectral criteria for the Parkes sample (Snellen et 
al. 2002), while the CORALZ sample (Snellen et al. 2004) 
was selected on the basis of radio morphology, not of the 
spectra. 

In Fig. |3 we compare the shapes of the radio spectra 
of the sources in the different samples, normalized to the 
peak frequencies and to the peak flux densities. It is ap- 
parent, and quantified in Table El that the median athick 
of Bolton et al. (2004) sources is considerably flatter than 
for the other samples. To homogenize this sample to the 
others, we have dropped sources with athick > —0.5. 

5. Results 

We have fitted separately the redshift and the observed 
peak frequency, v p ,o, distributions of GPS galaxies in the 
samples described above. The best fit values of the pa- 



Table 10. Best fit values of the model parameters for GPS 
galaxies. 



log(no) [Mpc- 5 (dlogL,)- 1 GHz -1 ] 


-14.13 


P 


0.83 


log(^,i) (GHz) 


2.98 


V 


1.20 


A 


0.0125 



ramctcrs have been obtained minimizing the chi-square 
function with the Minuit package (CERN libraries). We 
need to determine: 

— the parameters characterizing the luminosity func- 
tion, i.e. the normalization uq and the slope /3; the 
minimum and maximum luminosities, L m j n and L max 
respectively, both referred to the frequency of 300 
MHz, below the conventional minimum peak fre- 
quency of GPS sources, have been set at the minimum 
and maximum observed values: L max (300MHz) = 5 • 
10 35 ergs -1 Hz" 1 , L min (300MHz) = 10 29 ergs -1 Hz -1 ; 

— the slope, rj, of the power-law dependence of the emit- 
ted radio power on source age; 

— the slope, A, of the power-law dependence of the peak 
frequency on source age; 

— the initial value of the peak frequency, v v ^\ 

— the redshift of formation of the first peaked spectrum 
sources, z/; 
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Fig. 4. Comparison of the redshift distributions yielded by the best fit model (dashed) with the observed ones (solid). 



— the parameter k characterizing the cosmological evo- 
lution of the luminosity function. 

We have checked that the data do not require cosmologi- 
cal evolution of the luminosity function of GPS galaxies, 
confirming the finding of Dc Zotti et al. (2000), and we 
have therefore set k = 0. We thus have L+(z) — const = 
L a = 10 32 ergs^ 1 Hz~* [Eq. ©]. The fit is also insensitive 
to the value of z/, provided that it is larger than the maxi- 
mum estimated redshift of GPS galaxies in the considered 
samples. We have therefore fixed Zf = 1.5. 

Figures 0] and [S] compare the model with the observed 
redshift and peak frequency distributions. The overall 
agreement is reasonably good, indicating that the under- 
lying scenario is consistent with the current data. 

The best fit values of the parameters are given in 
Table ^| The formal errors on them derived from the x 2 
statistics are rather small, but we regard them as unrealis- 
tic in view of the many uncertainties due to the difficulties 
in the sample selection; also the redshift distributions are 
largely built using photometric redshift estimates rather 
than spectroscopic measurements. Thus we prefer not to 
report uncertainties that are likely to be deceitful. That 
the real uncertainties are probably large is indicated by a 
comparison of the present best-fit values of the parameters 
with those found by De Zotti et al. (2000) analyzing some 
of the samples considered here. In particular, the values 
of 77 and A are widely different [the difference of no fol- 
lows from that of A, see Eq. (J5J]. Still both analyses find 
positive values of 77 and negative values of p = — ?y + aA, 
implying that both the emitted radio power and the peak 



luminosity decrease with increasing source age, at vari- 
ance with the evolution models by Snellen et al. (2000) 
and Alexander (2000). 

The slope of the luminosity function is more stable: we 
find (3 = 0.83 while De Zotti et al. (2000) found /3 = 0.75. 
Both values are not far from the slope of the low lumi- 
nosity portion of the luminosity function of steep spec- 
trum radio sources, f3 ss ~ 0.69 (Dunlop & Peacock 1990; 
Magliocchetti et al. 2002). 

6. Discussion and conclusions 

Although new samples have substantially improved the 
coverage of the luminosity-redshift-peak frequency space 
of GPS galaxies, the assessment of their evolutionary prop- 
erties is still difficult. First of all, there is no fully agreed 
set of criteria to ascertain whether a source is truly a GPS. 
The frequently adopted GPS identification conditions rely 
on the spectral shape. However, as mentioned in Sects. [21 
and 01 different spectral requirements have been adopted 
by different groups and recovering a homogeneous set of 
data is not easy and may be even impossible. Other prop- 
erties that GPS sources should have, include a compact 
structure (size ^ lkpc), low polarization, low variability, 
and sub-luminal component motions. Regrettably, mea- 
surements of these quantities are available only for a lim- 
ited number of GPS candidates. The samples we are using 
may therefore be contaminated by sources of different na- 
ture. 

Also, a significant fraction of objects are still uniden- 
tified, and spectroscopy of identified objects is highly in- 
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Fig. 5. Comparison of the peak frequency distributions yielded by the best fit model (dashed) with the observed ones 
(solid). 



complete, so that many redshifts are estimated from op- 
tical magnitudes. Although GPS galaxies seem to have 
a rather well defined redshift -magnitude relationship, its 
dispersion is significant and may increase with redshift. 

With these premises we cannot expect to be able to 
come out with a clear assessment of the evolutionary prop- 
erties of GPS galaxies. Still, some interesting conclusions 
can be drawn. First, the simple luminosity evolution sce- 
nario for individual sources outlined in Sect. appears to 
be fully consistent with the data. We note however that, 
although the formalism stems from the self-similar evolu- 
tion models by Fanti et al. (1995) and Begelman (1996, 
1999), according to which both the parameters r\ and A 
are determined by the slope n of the density profile of 
the ambient medium, the data can be satisfactorily fitted 
only if the two parameters are treated as independent. 
This is not surprising, since self-similarity can be easily 
broken under realistic conditions, and indeed deviations 
from self-similarity were found in the two-dimensional hy- 
drodynamical simulations of Carvalho & O'Dea (2002). 

The fit is obtained for a positive value of rj and a neg- 
ative value of p, implying a decrease of the emitted power 
and of the peak luminosity with source age or with de- 
creasing peak frequency, at variance with the Snellen et 
al. (2000) model. On the other hand, our analysis con- 
firms the rather flat slope of the luminosity function, found 
by Snellen et al. (2000) who also report indications of a 
high luminosity break, not required by the data sets we 
have used. Snellen et al. (2000) argue that, in the frame- 
work of a scenario whereby GPS sources increase their 



luminosity until they reach a size ~ 1 kpc and dim there- 
after, during the CSS and extended radio source phases, 
the luminosity function of GPS galaxies can evolve into 
that extended radio sources. Our results suggest that the 
GPS galaxies are the precursors of extended radio sources 
with luminosities below the break of the luminosity func- 
tion (ibrcak ~ lO^ergs -1 Hz -1 at 2.7 GHz, cf. Dunlop & 
Peacock 1990, De Zotti et al. 2005). Our best fit model 
implies that the source luminosity is, roughly, inversely 
proportional to the source age. If GPS sources have typi- 
cal ages of ~ 10 3 yr and extended radio sources of ~ 10 7 yr, 
we expect the latter to be ~ 10 4 times less powerful than 
the former; as shown by Fig. ^ the maximum value of 
the peak luminosity we found is - lO^ergs^Hz -1 (but 
this may be a lower limit, since we have arbitrarily set 
the maximum redshift at Zf = 1.5), i.e. about 3 orders 
of magnitude larger than Lbreak- The GPS sources in the 
samples considered here are thus expected to evolve into 
large-scale radio galaxies with L < I/break- 

It must be stressed, however, that the uncertainties are 
very large, so that firm conclusions must await for larger 
samples and more complete redshift information. In par- 
ticular, we have checked that still acceptable fits can be 
obtained setting the exponent 7/ of the relationship be- 
tween the emitted radio power and the source age to the 
value implied by Begelman's model for n = 2, i.e. r; = 0.5 
(the total x 2 increases by 8\ 2 — 2). In this case, the source 
luminosities would typically decrease by a factor of ~ 100 
from the GPS to the extended radio source phase, and 
GPS sources could be the progenitors of extended radio 
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sources with luminosities both above and below L break- 
It is plausible that sources have a distribution of values 
of 77 whereby the lower values become increasingly rare, 
to account for the steepening of the luminosity function 
above I/break- Alternatively, the fraction of prematurely 
fading sources may increase with luminosity. If r\ = 0.5, 
the best fit value of A is 1.02, and the peak luminosity, 
L v {v p ) increases, albeit slowly, with decreasing peak fre- 
quency [Lp{u p ) oc Vp 25 , cf. eq. ©]• The other values of 
the parameters keep values very close to those in Table ITTTl 

Another key observable quantity potentially providing 
crucial constraints on evolutionary models is the linear 
size. The self-similar evolution models predict its depen- 
dence on source age, given the density profile of the ambi- 
ent medium, and on the other fundamental observable, the 
radio power. However, as noted above, we could not fit the 
data with strictly self-similar models, so that we no longer 
have a well defined relationship between radio power and 
linear size. Such a relationship, however, can be recov- 
ered through the observed relationship between turnover 
frequency and linear size, v p oc l~ s , with 5 ~ 0.65; we 
have L oc l~ n5 / x . The anticorrelation between the latter 
quantities implies that, since the radio power decreases 
with decreasing peak frequency, it decreases with increas- 
ing linear size. Once again, however, the uncertainties are 
very large. The formal best fit values of the parameters 
would imply an unrealistically steep decrease of the radio 
power, L, with increasing linear size, I, but, if r; = 0.5, 
L oc l~ - 32 and L p (v p ) oc P s / X oc Z 017 . 

In turn, fixing the exponent A of the relationship be- 
tween peak frequency and source age to the value implied 
by Begelman's model for n = 2, i.e. A = S = 0.65, we get 
again an increase of the minimum \ 2 by Sx 2 — 2, the best 
fit value of 77 is 0.76, while the other parameters are es- 
sentially unchanged. In this case we have Lp{v p ) oc ^°' 42 , 
Locr a76 ,andL p (i/ p )ocr - 27 . 

These examples illustrate the importance of the de- 
termination of linear sizes for large complete samples of 
GPS sources. Given their extreme compactness and their 
redshift distribution, images with milli-arcsec resolution 
are generally required. It is thus not surprising the the 
currently available information is scanty and inhomoge- 
neous. New VLBA data are however being acquired and 
analyzed, for example on the HFP sample (Orienti et 
al. 2005). Future evolutionary models should comply also 
with the distributions of linear sizes and with relationships 
of sizes with radio power. 

Complementary information on the proposed evolu- 
tionary sequence linking GPS sources to large scale ra- 
dio sources is provided by CSS sources, which could cor- 
respond to the intermediate phase when sources have 
expanded out of the narrow line region but are still 
within the host galaxy. More complete evolutionary stud- 
ies should take into account also these sources. 
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